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Abstract 



Oh! 

Gluon fragmentation to x CJ (lP) followed by single photon emission represents the 



dominant source of prompt J/^'s at the Tevatron for p±>6 GeV. Since fragmenting 
gluons are approximately transverse, their products are significantly polarized. We find 
that gluon fragmentation populates the helicity levels of Xdi Xc2 and J/^ according to 
D ( fc=o) : D (im=i) ~ 1 : 1, D (h =o) : D { \ h \=i) : D m =2) ~ 1 : 2.9 : 6.0 and D J/xs , ih=0) : 
Dji^{\h\=i) — 1 : 3.4. We also speculate that gluon fragmentation to the radially excited 
X C 2(2P) state followed by subsequent radiative decay could represent a large source of 
t(^'(2S) , s and potentially resolve the ip' deficit problem. A measurement of these states' 
polarizations would test this idea. 
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The production of the J/ty charmonium bound state is currently under active study 
at Fermilab 0. Until recently, the dominant sources of J/^'s at a hadron collider were 
believed to be parton fusion and B meson decay. These two processes respectively produce 
prompt and delayed J/W$ which can be distinguished via B meson vertex displacement 
measurements. Comparison between the theoretical prediction and experimental mea- 
surement of the transverse momentum differential cross section da(pp —>■ J/^ + X)/dp± 
section reveals that parton fusion alone underestimates the prompt J/W production rate 
at high transverse momenta by approximately an order of magnitude Such a large 

discrepancy between theory and data clearly indicates that another prompt production 
mechanism must be at work. 

Within the past few years, parton fragmentation has been examined as an alternate 
source of J/^'s Although fragmentation takes place at higher order in perturbative 
QCD than quark or gluon fusion, the falloff of the former with increasing transverse mo- 
mentum is much slower than that of the latter. So for p±>6 GeV, parton fragmentation 
represents the dominant source of prompt J/^'s. 

The first charmonium fragmentation functions to be calculated were D g ^j/q,(z) and 
D c ^>J/y( z ) which specify the probability for gluons and charm quarks to hadronize into 
J/ty as a function of its longitudinal momentum fraction z The only nonperturba- 

tive piece of information which enters into the lowest order computation of these S-wave 
fragmentation functions is the square of the charmonium bound state's wavefunction at 
the origin. The remainder of the calculation is based upon perturbative QCD. More re- 
cently, the fragmentation functions for gluons and charm quarks to hadronize into the 
lowest lying P-wave charmonium bound states x c o, Xci an d Xc2 have been computed 0. 
These Xcj states radiatively decay down to J with the branching ratios 0.7%, 27% and 
14% for J = 0, 1 and 2 respectively. After folding together these branching ratios with 
the P-wave fragmentation functions, one finds that gluon fragmentation to Xa followed 
by single photon emission to J/\& dominates at high p± over all other prompt mechanisms 
by more than an order of magnitude. When this J/\& source is included, the theoretical 
prediction for da(pp — > J/W + X)/dp±_ at y/s = 1.8 TeV moves to within a factor of two 
of recent CDF data. 

Most of the fragmentation functions which have been calculated to date describe the 
production of unpolarized quarkonium. However, it is straightforward to compute polar- 
ized fragmentation functions as well. Charm fragmentation into transverse and longitudi- 
nal J/^'s was considered in refs. and || and found to yield essentially no polarization. 
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J/vI/'s produced at a lepton collider like LEP are therefore not expected to be significantly 
polarized. But since gluon fragmentation to Xcj represents the dominant source of J/^'s 
at a hadron machine like the Tevatron, the polarized fragmentation functions D g ^ x (\h\)(z) 
where \h\ < J denotes the helicity of the produced Xcj must be determined before the de- 
gree of J/W polarization can be estimated. We present the results for these fragmentation 
functions and the Xcj and J/^> polarizations which they induce in this letter. 

To begin, we adopt the notation and general methods for computing P-wave frag- 
mentation functions established in 0. The lowest order diagrams that contribute to Xcj 
fragmentation are illustrated in fig. |I[ They may be evaluated using standard Feynman 
rules for quarkonium processes 0. The kinematic regime in which these graphs become 
important occurs when the lab frame energy qo of the incoming off-shell gluon g* is large, 
but its squared four-momentum q 2 is close to the square of the charmonium bound state's 
mass. We therefore neglect terms which are subdominant in the ratio q 2 /q 2 - The sum of 
the two diagrams in fig. [I] logarithmically diverges in the limit z — > 1 when the bound state 
carries off all the original gluon's energy and the outgoing gluon undergoes zero recoil. 
This infrared divergence is canceled by the diagram in fig. |2| which depicts the conversion 
of g* into a color-octet 3 Si state. The colored state can subsequently emit a soft gluon 
and decay to a color-singlet Xcj- The graphs in figs. 1 and 2 must be added together to 
obtain an infrared finite result. 

To determine the polarized fragmentation functions D uh\)(z), we need the polar- 
ization sums for the individual helicity levels of Xcj- These spin sums may be conveniently 
expressed in terms of an auxiliary light-like vector n = (1, —p) where p denotes a unit 
vector oriented along the three-momentum of the Xcj in the lab frame. The longitudinal 
and transverse polarization sums for the spin-1 boson can then be simply written in the 
covariant forms |10[ 

h=0 



a) 



\h\=i 



where 



1 p 2 

*aB — ~9a(3 H 5— 

p A 



(2) 



represent two-dimensional transverse and three-dimensional projection operators. The 
corresponding spin-2 polarization sums are given by 
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After a straightforward computation, we obtain the polarized x CJ fragmentation 
functions: 
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1 We used the high energy physics Mathematica package Feyncalc to perform much of the 
tedious algebra [11]. 
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In these expressions, M ~ 2m c denotes the charmonium bound state's mass, A represents 
the infrared cutoff, and H 1 = 72\R[(0) \ 2 /(nM 4 ) and H' 8 (A) = 8\R ( Q 8) (0)| 2 /(3ttM 2 ) re- 
spectively contain the squares of the derivative of the color-singlet P-wave and color-octet 
S'-wave radial wavefunctions. If the polarized fragmentation functions are summed over 
their helicity levels, we recover the unpolarized P-wave fragmentation functions reported 
in ref. 0. 1 

We adopt the parameter values M = 3500 MeV and a s {M) = 0.24. Following 0, we 
also set the infrared cutoff to A = M/2 and take H 1 ~ 15 MeV and H' 8 (M/2) ~ 3 MeV. 
The functions in eqns. (|]) - (|6|) can be evolved from the charmonium scale to higher 
energies using the Altarelli-Parisi equation and folded together with the gluon cross section 
da(pp — > g + X)/dp± to obtain the transverse momentum distribution of Xcj's produced 
at the Tevatron: 
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Since the gluon cross section is a very rapidly decreasing function of p±, it is a good 
approximation to evaluate the integral in (|7|) retaining just the terms proportional to 



The longitudinal and transverse Xci fragmentation functions were calculated in ref. [12]. Our 



results in eqns. (5a) and ( pq ) agree with those in [12] for z / 1 but differ for z = 1. 
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Longitudinally and transversely polarized Xci's are therefore produced at the Tevatron in 
the ratio D (h=o) : D qh\=i) — 1 : 1, while the helicity levels of Xc2 are populated according 
to D ( h =o) : -D (|h|=i) : D 01=2) — 1 : 2.9 : 6.0. Comparing these J = 1 and J = 2 

X c 2 X c 2 X c 2 

ratios to their unpolarized counterparts D o=o) : -D (|h|=i) =1:2 and D (h=o) : -D (|h|=i) : 

^cl -^cl ^c2 Xc2 

-D (|h|=2) = 1 : 2 : 2, we clearly see that the Xci an d Xc2 states produced as a result of 

Xc2 

gluon fragmentation are significantly polarized. 

The source of this sizable Xcj polarization can be traced to the fragmenting gluon. 
If the gluon were on-shell, its polarization would be completely transverse. The extent 
to which g* is actually off-shell modifies this result by only 0(q 2 /q^) terms. This effect 
can be seen most simply in the diagram of fig. |2]. The color-octet 3 S± cc in the \ccg) 
Fock state must inherit the gluon's polarization in order to conserve angular momentum. 
The subsequent transformation of this state into Xci an d Xc2 populates their respective 
\h\ = 0,1 and \h\ = 0,1,2 helicity levels in the ratios 1:1 and 1:3:6. As the color-octet 
terms in eqns. (||) - (^|) numerically dominate over the color-singlet terms, this explanation 
accounts in large part for the Xa polarization which we have found. 

We now turn to consider the radiative decay Xcj J/^+j. Since the electromagnetic 
branching ratio for XcO is more than an order of magnitude smaller than those for Xci 
and x C 2, we neglect its contribution to J/W production. The invariant amplitudes for 
the remaining J PC = (1,2) ++ P-wave charmonium states to decay via El transitions to 
the J PC = (1) S'-wave state must be parity even, charge conjugation symmetric and 
electromagnetic gauge invariant. They can be written down by inspection: 



Altarelli-Parisi evolution approximately cancels in the ratios of fragmentation functions. We 
therefore neglect it in our polarization analysis. 



iA(xd(p) -> J/*(p — k) + 7 (fc)) = g^^k^^e^e 
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Given these amplitudes, we can determine the photon's angular distribution in the 
Xcj rest frame. Letting 9 denote the angle between the photon's three-momentum in this 
frame and the Xcj's three- momentum in the lab frame, we form the dimensionless ratio 

J D (M r(y v 

R (J \cos9) = ] 

h= 
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(xi7 - J/* + 7) 
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which is a convenient measure of the angular distribution of photons from polarized x CJ 's. 
The explicit dependence of upon cos 9 is given by 
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If Xci and Xc2 were unpolarized, then p = 2/3, a = r = 2/5 and i?^ J ^ would become 
independent of cos#. But the fragmentation results in eqn. @ imply p ~ 0.49, a ~ 0.29 
and r ~ 0.61 and yield 

l? (1) (cos#) ~ 0.47 [1 + 0.21 cos 2 9] (12a) 
i? (2) (cos 9) ~ 0.46 [1 + 0.30 cos 2 0] . (126) 

In principle, measurements of and i?*- 2 ** would determine the polarizations of Xci and 
Xc2- But in practice, it will be much easier to observe the average angular distribution 

2 

^D g ^ XcJ Br( XcJ -> J/* + 7 )i^)(cos0) 
# (avg) (cos 0) = — ~ 0.47 [1 + 0.25 cos 2 9] (13) 
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and extract an average Xcj polarization. 

The amplitude expressions in (^) can also be used to derive the polarization of the 
J/ty which is induced by its Xcj progenitor. The feeddown from the separate Xcj helicity 
modes to those of the J/^> is given by 

1 
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After inserting the radiative branching ratios and Xcj fragmentation probabilities, we find 
that longitudinal and transverse J/^'s are produced in the ratio Dj^( h =o) : f j/^(w=i) — 
1 : 3.4. Equivalently, the ratio of transversely polarized to total J/^'s equals ( ~ 0.77. 
This ratio may be measured in the leptonic decay J/ty — ► If represents the angle 

between the lepton three-momentum in the J/W rest frame and the three-momentum of 
the J/ty in the lab frame, then 



dF 



ty-+tT) 3 

a cos 6 _ 

r(^ -> £+£-) 4 



(^1 — ^) — (^1 — -c) cos 2 O ~ 0.46 [1 + 0.25 cos 2 6]. (15) 



J/ty production via gluon fragmentation to Xcj consequently induces a 25% shift in the 
lepton pair angular distribution relative to unpolarized J/^'s. 

The Xcj and J/^> polarizations that we have found represent model independent pre- 
dictions of QCD which can be experimentally tested. Verification of these results would 
provide nontrivial checks of the entire fragmentation picture of quarkonium production at 
large p±. 

We have so far considered the production of only the lowest lying n = 1 radial level 
charmonium states. However, fragmentation ideas can be simply applied to higher radial 
levels as well. In particular, gluon and charm fragmentation to i^'(2S) have been studied in 
ref. 0. Even when fragmentation is included along with direct production, the theoretical 
prediction for ip' production underestimates the number of ^"s which have been observed 
at the Tevatron by roughly a factor of 50. This large gap between theory and data strongly 
suggests that some important ip' production mechanism still remains to be included. 

It is important to recall that ip' is the heaviest cc bound state which lies below the 
DD threshold. Therefore, n = 1 Xcj states cannot radiatively decay to ip' ', but their 
n = 2 counterparts can. None of these latter states which lie above the DD threshold have 
been observed. Estimates for their masses yield M(x c o(2P)) = 3920 MeV, M(x c i(2P)) = 
3950 MeV and M(x c2 (2P)) = 3980 MeV @. These mass values kinematically allow the 
5-wave transitions Xco(2-P) — ► DD and % c i(2P) — > D*D to occur. We therefore expect 
the J = and J = 1 excited x CJ 's to be broad and to have negligible branching fractions 
to lower cc bound states. However, angular momentum and parity considerations require 
the analogous decays x C 2(2P) — > DD and % C 2(2P) — > D*D for the J = 2 state to proceed 
via L = 2 partial waves. Although we cannot readily compute by how much these P-wave 
decays will be suppressed, it is possible that the branching fractions for Xci (2-P) transitions 



7 



to charmonium states below DD threshold could be significant. @ If so, the experimentally 
measured branching fraction Br(xb 2 (2-P) — > T(2S) + 7) = 19% in the bb sector suggests 
that the corresponding fraction Br(x C 2(2P) — > %p'{2S) + 7) could lie in the few percent 
range. We estimate that a 5% branching fraction would enhance the theoretical prediction 
for ip' production by more than an order of magnitude. Such an enhancement would help 
resolve the ip' deficit problem. 

Our proposal is admittedly speculative. If this idea is correct, then a ip' should be 
accompanied by a photon resulting from \ci (2-P) decay. Since the Xc2 (2P) is polarized to 
approximately the same extent as its n = 1 counterpart, the photon will be distributed in 
angle according to R^(cos9) as specified in eqn. ( |126|) . Moreover, the induced polarization 
for ip' should be slightly enhanced relative to that of J/W since Xci(2P) does not feed down 
along with x C 2(2P). A measurement of these radially excited states' polarizations would 
therefore provide a test of this possible ijj' production mechanism. 



4 The 1 D2 charmonium state is forbidden from decaying to DD by parity. Moreover, its mass 
is predicted to lie below the D* D threshold [13]. This state is therefore narrow. However, its 
production is suppressed, and it is expected to have a very small branching ratio to ip'"/. The 
contribution of the 1 Z?2 state to ip 1 production is consequently negligible. 
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Figure Captions 

Fig. 1. Lowest order Feynman diagrams which mediate gluon fragmentation to color- 
singlet P-wave charmonium bound states. 

Fig. 2. Lowest order Feynman diagram which mediates gluon fragmentation to a color- 
octet S-wave charmonium state. 
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This figure "figl-l.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9408352vl 
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